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The Hamiltonian conservative system of two interacting particles has been considered both in clas- 
sical and quantum description. The quantum model has been realized using a symmetrized two- 
particle basis reordered in the unperturbed energy. Main attention is paid to the structure of chaotic 
eigenfunctions (EF) and to the local spectral density of states (LDOS). A remarkable correspon- 
dence has been found for the shapes of EF and LDOS in the energy representation, to their classical 
counterparts. Comparison with the Band Random Matrix theory predictions has revealed quite 
significant differences which are due to dynamical nature of the model. On the other hand, a partial 
agreement is found by inserting randomness " ad hoc " in the dynamical model for two-body matrix 
elements. This shows that, at least for small number of particles, care must be taken when classical 
correlations are neglected. The question of quantum localization in the energy space is discussed 
both for dynamical and random model. 



PACS numbers; 05.45+b 

I. INTRODUCTION 

Quantization of classically chaotic systems has been 
addressed, from the very beginning, both to conserva- 
tive and time-dependent systems. In the latter case 
the important phenomenon of dynamical localization was 
discovered, connecting a classical quantity, the diffusion 
rate, to the quantum localization length of the corre- 
spondent equilibrium distribution^. Instead, in the case 
of conservative systems, important steps have been done 
in establishing some distinctive features which mark a 
quantum chaotic system from a integrable one: let us 
mention for instance, the non-Wigner-Dyson statistics of 
neighboring level spacings, or the scarring of eigenfunc- 
tions along some classical periodic orbits. However, the 
possibility of quantum localization effects in such systems 
has been scarcely explored until recently, when a clue in 
this direction was foundo by investigating a particular 
class of random models: the Wigner Banded Random 
Matrices (WBRM) ensemble. 

For such an ensemble, whose introduction dates back 
to Wigner himselm, it is possible to obtain a series of re- 
sults which allow for a definition of quantum localization 
within the classical energy surface. These results, when 
extended to Hamiltonian systems, would impose severe 
quantum limitations on the behavior of classical ergodic 
systems. The important result is that room is left for 
quantum localization and this can obtained directly from 
the knowledge of local spectral density of states (LDOS) 
and eigenfunctions (EF). Quite surprising, both quanti- 
ties have well defined classical limits (see Ref.B), which, 
generally speaking, have received scarce attention before 
now. 

On the other hand, in order to acquire physical rele- 



vance, it is clear that such results should be extended to 
real physical systems, where the origin of randomness is 
purely dynamical. This we do, in this paper, by consider- 
ing a classically chaotic two-interacting spin system with 
a finite Hilbert space. Our purpose is analyzing the struc- 
ture of eigenstates and of LDOS, and comparing it with 
expectations based on previous Random Matrix studies 
and with their classical counterparts. 

First of all we find that, when written in the eigenba- 
sis of two non-interacting particles, reordered in the un- 
perturbed energy, the Hamiltonian matrix has an overall 
banded structure. About the shape of eigenfunctions and 
the LDOS, we find that our quantum results, in the av- 
erage, follow the behavior of similar quantities computed 
from WBRM only approximately at best. Nevertheless, 
they follow remarkably well the behavior of their classical 
analogs, that we actually compute in the present paper. 

On the other hand, the correspondence with Random 
Matrix Theories (RMT) is restored on artificially ran- 
domizing our Hamiltonian. The lesson we draw from 
this result is that, although RMT quite well reproduces 
fluctuation properties of spectra of real chaotic Hamil- 
tonians, some correlations are missing in their structure, 
which are essential in giving the correct semiclassical be- 
havior when detailed questions about the structure of 
eigenfunctions are asked. It is of course possible that 
a better correspondence with RMT will be restored with 
systems with a larger number of particles; for the time be- 
ing, however, our results appear to indicate that caution 
is needed in carrying over results from RMT to Hamilto- 
nians which have a smooth, well defined classical limit. 
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II. THE MODEL 



(5) 



pThe model has been proposed and widely investigated 
inP. Here we review few fundamental facts about its clas- 
sical and quantum behavior. It describes two coupled 
rotators, with angular momentum L and M with the fol- 
lowing Hamiltonian: 



H = A{L, 



(1) 



It may be used to describe the interaction of quasi- 
spins in nuclear physics or pseudo-spins in solid state 
systems. Choosing as the unit of time and AB~^ 
as the unit of angular momentum, it can be written as 
H = Hq + V where Hq = L,^ + and V = L^M^. The 
constants of motion arc H = E, and Af^. 



Fixing the values of L and M it can also be show: 



J 



that the total energy must be bounded: 

= (L2 + 1)(m2 + 1) 



(2) 



for LM >1. 

It is worth to mention that in this form the dynami- 
cal variables L, M are not canonical. On the other side, 
the usual Hamiltonian form, with the canonical variables 
QijPi, i — 1,2 can be recovered by means of the following 
transformation : 



Lx = \ JE^ -pI cosgi 
Ly = ^IP- -p\ sin gi 
Ez = pi 

= y^^P-pi cos 52 

My = ^M'^-pI sings 

= P2 



(3) 



keeping and as constantsQ. In these variables 
the Hamiltonian reads: 



H 



Pi + P2 + \J E^ - pi\J M'^ - pI cos qi cos q2 (4) 



The analysis of the surfaces of section reveals a large 
number of regular trajectories covering invariaiit tori 
when L'^ , are both very small or very largeQ. To 
simplify the problem we set L = M. In such a case 
the most interesting situation occurs when 1 < L < 10 
where, depending on the energy value regular and 
chaotic regions coexist. Typically when \E\ is close to 
Emax = L'^ + 1 trajectories are regular while for E ~ 
islands of stability become very small and chaotic motion 
dominates. 

Quantization follows standard rules, and angular mo- 
menta are quantized according to the relations — 
= n^l{l + l) where I is an integer number. Therefore, 
for given I the Hamiltonian is a finite matrix, and the 
semiclassical limit is recovered in the limit I — > oo and 
h ^ keeping constant. 

In our approach the Hamiltonian is represented in the 
two-particles basis \lz,mz) where the matrix elements 
have the form. 



and 

^/{l + -U + l){m + m,){m - + 1) 
with lz,mz integers, —I < lz,mz < I. 



(6) 
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FIG. 1. Structure of the Hamiltonian matrix in the sym- 
metrized basis for I = m — 11, (TV — 144) and Jz even. 
Dots indicate off-diagonal elements different from zero. The 
bandwidth b is maximal at the center where b = 21 + 1. Few 
different Ho shells are shown by the horizontal lines. 

The z-component of the total angular momentum 
Jz = Lz + Mz (which is the same as the unperturbed 
Hamiltonian Hq) obeys the selection rules AJz — 0, ±2?i, 
so the subspace spanned by the states with odd Jz/?i 
can be separated from that with Jz/h even (there are 
no matrix elements for the transition between them) . In 
what follows, we fix Jz/h = Hq/Ti even. As a result, 
the matrix describing the Hamiltonian has a dimension 
N — 2^^ -I- 2/ -I- 1. We have also to take into account 
the symmetry degeneracy with respect to the exchange 
of particles. Below, we consider only symmetric states. 

Let us now explain how the Hamiltonian matrix is con- 
structed. Once I is fixed, there are 21 + 1 single-particle 
levels (-/|, {-I + 1|, . . . , (-11, (0|, (1|, ...,(/ The 
ground state is represented by two particles in the low- 
est single-particle level, which we label as (— Z|; it 
has an unperturbed energy Eq — —2lh. The first ex- 
cited state is doubly degenerate and the two eigen- 
states having the same energy Ei = {—21 + 2)h are 
((-/,-/ + 2| + (-/ + 2,-/|)/V2 and {-l + l,-l + l\. The 
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former state corresponds to one particle in the single- 
particle level (— / + 2| and the other in the single-particle 
ground state The latter state, to two particles in 

the single-particle level (— / + 1|. We call i?o-shell the set 
of states having the same value of unperturbed energy. 

It is easy to prove that in the symmetrized basis each 
shell with Hq fixed and even has a degeneracy p — Z -I- 1 — 
\HQ\/2h, and the dimension of the Hamiltonian matrix is 
N ^ (l + 1)^ due to the relation 

Ho/2=-ni 

We then reorder the Hamiltonian matrix according to 
the increasing unperturbed energies and we call |n) the 
resulting two-particles symmetrized ordered basis. As 
a result, the off-diagonal matrix elements are symmetric 
with respect to the two main diagonals i?„^„ — (n|_ff|n) 

and HN-n,N-n- 

Diagonal matrix elements are constructed from the un- 
perturbed Hamiltonian Hq; they are given by the eigen- 
values —2lh, {—21 + 2)?i, . . . , 2lh and are disposed along 
the principal diagonal starting from the lowest left cor- 
ner. One should note that diagonal elements of the per- 
turbation V vanish due to (||). The global structure of 
the matrix -ff„i,n is shown in Fig.^ The next (to the 
principal one) diagonals Hn^n±i correspond to transi- 
tions inside each i?o-shell while the "arcs" connecting 
the two corners represent transitions between neighbor- 
ing shells having A_ffo = ±2?i. Such a global structure of 
the Hamiltonian matrix is not a peculiarity of this model 
but it corresponds to the so-called "shell model" reprei- 
sentation widely uspjcLin atomic and nuclear physicsB'Q. 
It was shown in Ref.l3illJ that, generic properties of eigen- 
functions in this basis, can be directly related to single- 
particle operators, in particular, with the distribution of 
occupation numbers for single-particle states. 

The Hamiltonian matrix has a clear band structure, 
with the bandwidth b ranging from 1 at the corners up to 
6 = 2/ + 1 in the middle. However, this structure differs 
strongly from that of staridard Wigner Band Random 
Matrices (see for examplea^El and references therein). 
Moreover, non-zero off-diagonal matrix elements are pos- 
itive and the mean and variance of the distribution of 
these matrix elements depend on the classical parameter 
~ h'^l{l + l) only. To be more precise, if one assumes a 
continuous distribution of the matrix elements, it can be 
shown (see Appendix 1) that = (v^) — (v)"^ ~ (i/4)^. 

There are of course semiclassical corrections to this 
estimate, but the variation of h in one order of magni- 
tude (available in our numerical study) change the ratio 
CT^/(L/4)^ by less than 1 %. In the same way the mean 
value can be semiclassically estimated as (w) ~ (ttX/S)^. 
The agreement between these simple semiclassical formu- 
las and our numerical data is shown in Appendix 1. 
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FIG. 2. Density of states for the total Hamiltonian Ho + V, 

(a) , and for the unperturbed one Ho, (b) as a function of the 
rescaled energy e = E/Emax for I = 39 and L = 3.5. For 
comparison, the fitting Gaussian (a) and the triangular curve 

(b) . 

The model is highly non-perturbative since the pertur- 
bation spreads the levels of the inner _ffo-shell all over 
the allowed energy range, see Appendix 2. As a re- 
sult, the perturbed spectrum is broader than the un- 
perturbed one; this is an effect of the non-zero mean 
value of off-diagonal matrix elements. While the unper- 
turbed spectrum has degenerate energy levels with spac- 
ing h (therefore, with density of levels po 1/h and 
spectral radius i?° = 2M ~ 2L), the perturbed one has 
(Z + 1)^ non-degenerate states within a spectrum of radius 
Ra ^ L'^ + I, which gives the density p ~ ~ Po/^- 

Numerical data show that the density of states changes 
from a triangular shape for the unperturbed Hamilto- 
nian (one more state is added/subtracted at each neigh- 
bor level for Hq negative/positive) to the Gaussian form 
which is known to be generic for realisitie finite systems 
like atoms and nuclei, see for exampleEll3. In Fig.^ the 
perturbed and the unperturbed density are shown for 
a typical case, together with the corresponding fitting 
curves. 



III. QUANTUM-CLASSICAL 
CORRESPONDENCE FOR EE AND LDOS 

The subject of this section is the analysis of the LDOS 
(also known as "strength hmction" or "Green spectra") 
and of the structure of eigenfunctions, together with their 
classical analogs. In the quantum description, all infor- 
mation is contained in the matrix constructed from the 
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eigenfunctions V-'n(£'m) of the total Hamiltonian H repre- 
sented in the ordered unperturbed two-particle basis |n). 
Here ipn{Em) is the n-th component of the eigenfunction 
having Em as eigenvalue. This matrix is assumed to be 
reordered in eigenenergies Em- 

In the classical limit the unperturbed energy Eq is 
not constant when the (chaotic) trajectory of the to- 
tal Hamiltonian H fills the H — E — const sur- 
face. Indeed it fills a range of values which are dis- 
tributed according to the ergodic measure on the con- 
stant enetgy surface, yielding a distribution function 
VF(i?o|-E)H. This distribution can be easily numeri- 
cally calculated taking a sample of chaotic trajectories 
u{t) = {L,{t),Ly{t),L,{t),M,{t),My{t),M,{t)) having 
the same fixed values of -E and = M^. Following these 
trajectories, one can calculate HQ{u{t)) = Lz{t) + Mz{t) 
taken at equal instants of time and find the distribution 
of Ho over the energy bandll3 defined by sojourn times. 

The quantum analog of this distribution is provided by 
the relation 

Wn{E) = {\MEm)\^)m (7) 

where the average (...) is taken over those eigenfunc- 
tions which have an eigenvalue Em in a fixed small energy 
interval around a given energy E. Such an average has 
been done in order to smoothen the fluctuations which 
affect individual eigenstates; we would like to note that 
for our dynamical model, unlike random matrix ensem- 
bles, there is no possibility of ensemble averaging. The 
distribution Wn {E) gives the average shape of eigenstates 
represented in the unperturbed two-particles basis |n). 

In order to obtain the quantum distribution W{Eo\E) 
one needs to switch to the unperturbed energy repre- 
sentation, n — !■ E'^. Technically this can be realized by 
introducing small energy bins AE and counting the cor- 
respondent probability within them, 

W{Eo\E) = J2 Wn{E)d{Eo - E°) (8) 

n 

Similarly, we can define the distribution w{E\Eo)- In 
the quantum case this distribution is the LDOS defined 
by 

w{E\Eo) = Y,{\MEm)\^)nS{E - Em) (9) 

m 

where the average is now taken over a number of val- 
ues of n, such that the eigenvalues E'^ belong to a small 
interval around the given unperturbed energy Eq. The 
presence of degeneracy in this case provides an obvious 
way of taking average. The corresponding classical func- 
tion can be found by noticing that the trajectory does 
not fill the whole surface Hq = Eq but is restricted to 
an invariant manifold specified by the value of m. Giv- 
ing equal weight to all m values corresponding to a given 
value of Hq , exactly matches the quantum averaging used 
in (0). Then the classical distribution can be evaluated 



analytically, since the classical unperturbed Hamiltonian 
Hq is integrable. Indeed, the unperturbed solution uo(t) 
for = and Hq = is given explicitly by : 



Kit) = 




cos 2(t — 


<P) 


Llit) = 




sm2(t — 


0) 


Lb) = 


m 






MSit) = 




m? cos 2(i — 




Mm = 
Mht) = 


VL^- 
— m 


m? sin 2{t — 






(10) 



where < < 27r, < ^ < 27r and \m\ < L depend on 
the initial conditions. Therefore, 

H{uQ(t)) = (L^ - r7i^)cos2(i- 0) cos2(i-^) 

This means that the classical distribution of H is given by 
Pbiu) where y — L'^{1 — xl) cos7ra;2 costtxs is a function 
of the random variables —1 < Xi < 1, i ~ 1,2,3. 
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FIG. 3. Typical eigenfunctions for the case L = 3.5, I = 39. 
a) The second excited state tpri{E2) (corresponding to a clas- 
sically integrable region) . b) Eigenfunction for the energy Em 
close to zero (middle of the spectrum, corresponding to the 
chaotic region). 

IV. STRUCTURE OF EIGENFUNCTIONS 

The quantum model has been already studied in Ref.i, 
but previous studies have not addressed the structure 
of eigenfunctions in the two-body particle basis. This 
representation is quite natural and corresponds to a well 
known procedure in the physics of interacting particles. 

In our dynamical model the structure of eigenfunctions 
strongly depends on their energy because in the classical 
limit for low and high energies (|e| = \E\/Emax ~ 1) the 
motion is regular while in the center of the energy band 
(|e| ~ 0) is chaotic. One can, therefore, expect that in 
the classical limit {h ~ L/l ^ 1) the eigenstates corre- 
sponding to regular or chaotic regions are very different. 
This is, indeed, clearly seen in Fig.^ where two eigen- 
states are plotted in the unperturbed two-particle basis 
for the second (from the bottom of the spectrum) eigen- 
statc (e ^ —1) and for the eigenstate chosen in the center 
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of the energy band (e ~ 0). Comparing these two (typi- 
cal) eigenstates, one can see that there are strong correla- 
tions between components of the "regular" one (Fig.^ a) 
while the "chaotic" eigenstate can be treated as random 
along the whole basis |n). We also would like to note that 
the regular eigenstate has many "principal components" 
and that it looks more or less extended. This again indi- 
cates that the perturbation is quite strong and effectively 
couples many unperturbed states. 

A much more accurate analysis of eigenfunctions is ob- 
tained by studying their localization lengths. Since the 
basis is finite and eigenstates can be extended along it, 
here we use different measures of localization lengths, 
basedi-op their entropy TL and participation ratio V see, 

.K.yM 



e.g 



Ih{E) = 2.08 exp{-7^} 



where 



kpr{E) = ■i/V 



N 



(11) 



W = 5] |V'„(i?)p;n|^„(i?)|^ 



n=l 



and 



N 



The normalizing coefficients 2.08 and 3 were chosen 



in order that Ih = Iv 



N in the limit case when all 



components ipniEm) are independent Gaussian random 
variables. Here N = (Z + 1)^ is the size of the two-particle 
basis. 

Further information can be extracted from the cen- 
troids ric of eigenstates and from their "widths" : 



nc{E)= Y.nAi^n{E)\' 

UE)^ {En\ME)\Hn-n,iE)rf' 



(12) 



In Fig.^ we present numerical results for the above quan- 
tities as functions of the rescaled energy e. 

First, we note that the entropy and inverse participa- 
tion ratio localization lengths, Ih and Upr, are approxi- 
mately equal and show the same behavior, namely, the 
delocalization along the whole basis in the middle of the 
spectrum, and the localization at the spectrum edges (see 
Fig.^ (a-b)). We also note that due to the underlying 
symmetry of the model the above quantities lipr, Ih are 
symmetric around e = 0. 

Even in the center of the spectrum, where eigenstates 
are in average maximally extended, there are big fluctu- 
ations in the value of localization lengths. This indicates 
that in the classically chaotic region there are some eigen- 
states which can not be treated as completely random 
and delocalized over the energy shell. A careful study 
shows that such eigenstates are characterized by an ex- 
tended background with some pronounced peaks (the so- 
called "sparse eigenstates"). Such eigenstates may result 



in the absence of equilibrium and iiVithe lack of standard 
statistical description, see details irJj. 



(a) 



(b) 
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£ £ 
FIG. 4. Measures of localization lengths for eigenfunctions 
versus the rescaled energy for the case L = 3.5, I = 39. a) 
entropy localization length Ih, b) localization length hpr de- 
fined as the inverse participation ratio, 
variance, la, d) centroid ric 



square root of the 



One can also see a clear regular structure in the de- 
pendence of Ih and kpr on the energy at the edges of the 
spectrum, which reflects the regular character of eigen- 
states. 

The other two quantities, ric, la give an information 
about the "position" and "width" of eigenfunctions in the 
two-particle basis, see Fig.||(c,d). In contrast to Ih and 
lipr the "width" 1^ reveals a quite unexpected minimum 
at the center of the spectrum. Additional numerical anal- 
ysis shows that this is a result of different "sparsity" of 
chaotic states depending on the energy. Namely, chaotic 
eigenstates are more compact at the center of the en- 
ergy band than far from it. In fact, the ratio la/lH,ipr 
can be used to extpajCt an information about the sparsity 
of eigenstates (seell3). Indeed, two eigenstates with the 
same value of lH,ipr can have very different values of la 
depending on whether principal components (those with 
relatively large values of ipniEm)) are clustering around 
some center (small la) or randomly scattered over the 
whole unperturbed basis (large la)- 

Additional information can be obtained from the de- 
pendence of the centroids of eigenstates on the energy, 
see Fig.^ (d). Apart from fluctuations and excluding the 
regular part, this dependence is linear, which means that, 
in average, centers of eigenstates are located at the center 
of the energy shell covered by the classical distribution 



5 



'w{E\Eo). This generic feature of chaotic eigenstatea has 
been studied in greater details in WBRM modelaJ. In 
particular, it was shown that those eigenstates, which are 
completely extended in the whole energy shell, are char- 
acterized by maximal statistical properties of quantum 
chaos. For example, in that case the statistics of energy 
spectrum follows the predictions of Random Matrix The- 
ory, such as the Wigner-Dyson form of the distribution of 
spacings between neighboring energy levels. On the other 
hand, localization of eigenstates within the leaergy shell 
leads to the so-called intermediate statisticalj (which is 
intermediate between the Wigner-Dyson and the Poisson 
statistics). 



not in the whole unperturbed basis !). This means that 
the observed scattering of the centroids of eigenstates 
(see Fig.||(d)) is, in fact, quite weak and does not lead to 
noticeable localization in the energy space. 

Interestingly, the size of these ergodic eigenfunctions 
is smaller than the total energy band. The distribution 
Wq is infact restricted between the minimum and the 
maximum value that the function Hq = Lz + can 
assume under the constraints E — Lz + Mz + L^M^, 
Ll + + Ll = + + = . It can be easily 
proved, using the Lagrange multipliers methods that for 
E = 0, \Eo\ < 2{\JlP- -I- 1-1), which in Fig.|| corresponds 
to \SuppW(tQ\e)\ < e~ 0.398 
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FIG. 5. a) Shape of eigenfunctions in the energy represen- 
tation (full line) and classical distribution M^(eo|e) (dashed) 
for the case L = 3.5, I = 39, obtained by averaging over the 
central / -I- 1 eigenfunctions with eo = 0; b) the same plot in 
semi-log scale. 

Such a localization is reflected in the fluctuations of 
Tic around the center of the shell (linear dependence on 
e). Indeed, if eigenstates are localized, their centers tIc 
are typically scattered within the energy shell leading to 
strong fluctuations of ric; instead, this cannot happen if 
they fill the whole energy shell. It is important to stress 
that localization in the energy shell is different from that 
in the unperturbed basis, as found from Eqs.([ll]). 

The result presented in Fig.^(d) shows that localiza- 
tion in the energy shell, if any, is quite weak. A more 
direct analysis of the degree of localization in the energy 
shell is provided by the direct comparison of the average 
shape of eigenstates in energy representation to its clas- 
sical analog. The results are presented in Fig.|^ where 
quantum and classical W{Eq\E) (see Eq|8|) are plotted 
versus the rescaled energy e = EjE^ax- 

One can see that the only important difference is a 
sort of a weak quantum tunneling in the classically for- 
bidden region (the tails of the classical distribution are 
sharper than the quantum ones). Anyway, the good cor- 
respondence between quantum and classical distributions 
shows that for the chosen parameters the model is in a 
deep semiclassical region and, globally, the eigenstates 
should be treated as ergodic ones (in the energy shell. 
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FIG. 6. a) The BS for an n value corresponding to a shell 
edge indicated as a dashed line in Fig. (8). b) The BS in the 
middle of //o = shell. 



V. LDOS STRUCTURE 

Of special interest is the structure of LDOS which is 
widely discussed in many applications in atomic, nuclear 
and solid states physics. The importance of this quantity 
relates to its physical meaning: it shows how an unper- 
turbed state |n) "decays" into other states due to in- 
teraction. In particular, the inverse width of the LDOS 
is associated with the mean "lifetime" of a chosen basis 
state. 

As was indicated above, the LDOS structure can be 
extracted from the matrix ?/>„(_£„) by fixing an unper- 
turbed state |n) and searching the dependence on m. 
Therefore, we can adopt the same procedure as we did 
above when analyzing the structure of eigenstates. In 
comparison with Fig.|3| we show two such "matrix lines" 
corresponding to basis states (BS), with close n values, 
taken from the center of matrix (see FigJ^). In fact, the 
BS-lines in this matrix correspond to the expansion of un- 
perturbed (basis) states in the exact eigenstates. While 
the structure represented in Fig.|^(b) is typical, the form 
of the BS in Fig.^(a) is only observed around some spe- 
cific value of n. 

To better understand the meaning of these peculiar 
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n— values we have computed the locahzation length (com- 
pare with Eq.(|l])) 

kprin)^ 3/Ern\ME,n)\' (13) 

in some range of n. The data presented in Fig.^ reveal 
a global periodic structure of BS, from which one un- 
derstands that the peculiarity of BS reflected in Fig.^ 
results from the degeneracy (inside each shell) of the un- 
perturbed spectrum. It is then convenient to consider 
in the following analysis, as a reference, the central shell 
Ho = only (set of Z-K 1 BS ). 
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FIG. 7. Dependence of the inverse participation ratio of the 
n— th BS of the H matrix as a function of n. Here is L = 3.5, 
/ = 39. The Ho — shell is inside the two vertical dashed 
lines. 

Now we discuss the structure of the BS in energy 
space which is, in fact, the LDOS (Eq.^). According to 
result^y obtained for WBRM, the shape of LDOS typ- 
ically changes from the Breit-Wigner (BW) law to the 
semicircle when an effective perturbation is increased. In 
particular, the BW is expected when 

1/V2^ < poV" < 7^2^ (14) 

where po is the density of the unperturbed spectrum and 
b is the effective band width of the Hamiltonian matrix. 
The first inequality is related to the non-perturbative 
character of the coupling (which is always verified in this 
model) while the last, rewritten as 

27rpoV^ = Ff < b/po, 

simply means that the spreading width Tp of such distri- 
bution has to be much smaller than the energy bandwidth 
b/ pQ. If we formally apply the above conditions to our 
case, we get (in units of e = E/Emax) the following re- 
lations: b/po ~ l/L and 2t:pqV^ ~ {-k I2){L/A)H. With 



our data the second condition in Eq.(|j) is strongly vi- 
olated. In fact the random matrix agrgument leading to 
the above results rests on the assumption that the band 
in the Hamiltonian matrix be "full" ; in our case, instead, 
we have a large sparsity (many vanishing matrix elements 
inside the band). 
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FIG. 8. a) Quantum LDOS distri- 

bution w(e|eo) = w{E / Emax\Eo/ Emax) (fuU line) and best 
fitted Breit-Wigner distribution in the range specified above 
(dashed line) for the ceise L = 3.5, I = 39, obtained averaging 
over / -I- 1 values of BS for the Ho = shell; b) the same plot 
in semi- log scale. 

In Figjs] we show the structure of the LDOS for the 
BS corresponding to the center of the unperturbed spec- 
trum, in comparison with the Breit-Wigner fit which is 
performed within the interval {—b/2po,b/2po). 
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FIG. 9. a) The LDOS distribution w;(e|eo) : quantum (full 
line) and classical (dashed line) for the case L = 3.5, / = 39, 
obtained by averaging over I + 1 values of BS for Ho = shell; 
b) the same plot in semi-log scale. 

One should stress that outside of the energy interval 
corresponding to the band sizfij— the tails of LDOS are 
known to be highly non-genericQiid'Ej'Ll depending on spe- 
cific properties of the model. As one can see, inside this 
energy interval the shape of LDOS can be roughly asso- 
ciated with the BW form. On the other hand, outside. 
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the tails decay very slow, compared to those given by the 
BW. Such a fown of the tails is also different from the case 
of the WBRMEl where outside the band energy range the 
tails decay extremely fast (even faster than exponential) . 
In general, the above results seem to indicate that the 
effective perturbation corresponds just to the condition 
when the BW approximation starts to fail. 

Let us now compare quantum and classical LDOS. In 
Fig.^ we give an example of such distributions in the 
whole energy shell. They coincide with a high accuracy, 
apart from the regions very close to the energy shell edges 
(where quantum tunneling is significant). 

This again means that the system is in a deep semi- 
classical regime. We remark that, in this model, going to 
the quantum regime in the chaotic energy region calls for 
very small matrices, for which fluctuations are extremely 
strong. 




in order to have the same energy range as for it;(e|eo): 
W^(eo|e) = iW{eo/e\e/e) where e = 2{^/IJTl- 1). The 
rescaled distribution Wr is presented in Figjl^ together 
with the distribution i/7(e|eo). After such a rescaling both 
distributions coincide quite well which again indicates the 
absence of the localization. 
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FIG. 11. Typical eigenstates for the case L — 3.5 
and / = 39 and Gaussian random non-zero off-diagonal el- 
ements, (a) and (b) have close eigenvalues in the middle of 
the spectrum, respectively e — E/Emax = 1-48 ■ 10~^ and 
e = E/E^ax = 2.295 • 10"^ but very different inverse partici- 
pation ratios ( kpr/N — 0.0027 for the first and hpr/N = 0.87 
for the last, where A'' is the matrix size.) 



VI. RANDOM TWO-BODY INTERACTION 



FIG. 10. a) Classical LDOS distribution w{e,eo) (full line 
) and the rescaled classical distribution VFr(eo|e) for the case 
L = 3.5, / = 39. b) the same plot in semilog scale. 

An important question is the relevance of the shape of 
eigenstates to that of the LDOS. As was shown in the 
model of WBRMLO in some range of parameters (for not 
very strong perturbation) the two shapes are very close to 
each other, which is a manifestation of the ergodic struc- 
ture of eigenstates in the energy shell. On the other hand, 
with an increase of perturbation the LDOS was found to 
tend to the semiciixle, for which strong localization turns 
out to be possibleo. This localization manifests itself in 
different average shapes of the EF and LDOS. Namely, 
the width of the EF in the energy representation is less 
that the width of the LDOS; the latter defines, in fact, 
the width of the whole energy shell. 

Direct comparison of Fig.(||) and Fig.(^) in our dy- 
namical model shows a remarkable different energy range 
for the LDOS and EF distribution. As was discussed 
above, the energy width of eigenfunctions in the semi- 
classical region (in the energy representation), is much 
smaller than the width of the spectrum because it is sub- 
ject to an additional constraint. We can take into ac- 
count this restriction and rescale the distribution VF(eo|e) 



In this Section we modify our model of two interacting 
particles by assuming a completely random interaction 
which preserves some global properties of the original 
dynamical model ([^). Namely, the unperturbed part Hq 
is taken exactly the same as in the dynamical model. 
However, we replace non zero matrix elements of the dy- 
namical model with random and independent variables. 
Moreover, we choose a Gaussian distribution of these ran- 
dom matrix elements with the same mean and variance 
as for the dynamical elements. In such a way we can 
reveal the influence of dynamical correlations which are 
due to the specific form of the interaction V. Below we 
follow the same procedure described in previous sections 
when studying eigenstates and LDOS. 

Numerical data for the "randomized" model show that 
global spectral properties are the same as in the dynam- 
ical model. Namely, the perturbed spectrum is enlarged 
with respect to the unperturbed one and the density 
of states keeps the same Gaussian shape with the same 
mean and variance. 

On the other hand, the analysis of eigenfunctions re- 
veals clear differences. Typical shapes of eigenstates and 
BS are shown in Figjl|. Comparing with the correspond- 
ing Fig.^(b), one notes that extended states look chaotic, 
similar to those found in the dynamical model. However, 
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differently from the dynamical model, a few strongly lo- 
calized states now appear even in the center of the energy 
band. A typical example of such an eigenstate is given 
in Fig.|l|(a). 

To analyze the global characteristics of all eigenstates 
we have calculated different localization lengths Ih and 
lipr as well as the width and the centroids ric accord- 
ing to Eq.(0|2|). The data reported in Fig.|l^ should be 
compared with those in Fig.^ for the dynamical model. 
As expected, for the random model there are no corre- 
lations in the energy dependence for large/small energy 
(|e| ~ 1), compared with Fig.^. However, close to the 
edges of the energy spectrum, the eigenstates can not be 
treated as chaotic since the number of "principal compo- 
nents" in such eigenstates is quite small (this is revealed 
by small values of localizations lengths Ih, hpr, la)- This 
is a result of the perturbative localization which typi- 
cally occurs for states close to the ground state. In what 
follows, we exclude such states from our consideration. 



structure of eigenstates are much stronger than in the 
dynamical model. 






FIG. 13. (a) The EF distribution for the Gaussian random 
case with L — 3.5, I — 39, obtained by averaging over / -I- 1 
central eigenfunctions; (b) the same as (a) in semilog scale. 

Despite these fluctuations, the global structure of the 
EF seems to remain the same. This is marked once more 
by the distribution of eigenfunctions in the energy space 
(the analog of Fig.^ is now shown in Fig.^Sj) which has 
the same shape as in the dynamical model. However, the 
LDOS distribution for the random model shows striking 
difference, see Fig.|lJ. Indeed it can be described, apart 
from the central peak, by the semicircle lawH'OH. 



(c) 
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FIG. 12. Measures of localization lengths for eigenfunc- 
tions for the Gaussian random case L = 3.5, I — 39; compare 
with Fig. (5). 

For chaotic eigenstates, the various measures of local- 
ization lengths give average values typically less than in 
the dynamical case. This holds especially in the middle 
of the spectrum _E ~ where we can now find a rela- 
tively large number of sharply localized eigenstates (with 
hpr/N ^0.1 where N is the matrix size), see Fig.O (a). 
In the same Fig.|l^(d) one can also observe a much more 
stronger scatter of the centroids of eigenstates transverse 
to the diagonal (compare with Fig.^(d) of the dynami- 
cal model). The same features have been found for ba- 
sis states. These data indicate that fluctuations in the 



FIG. 14. (a) The LDOS distribution for the Gaussian 
random case with L — 3.5, / — 39, obtained by averaging over 
I + 1 central lines; (b) the same as (a) in semilog scale. 

This surprising result is quite significant in the ligt 
application of Random Matrix models. Indeed, ir£f 
it was found that localization in energy shell for the 
WBRM may occur only when the LDOS is characterized 
by the semicircle law. Therefore, the important ques- 
tion is whether the semicircle law is a quantum feature 
or it can also occur in classical dynamical conservative 
systems. What we have found here is that the semicircle 
law has nothing to do with the semiclassical limit in our 
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model. It seems to be dictated by quantum randomness 
rather than by the pseudo-randomness resulting from the 
classical chaos. 

By comparing the shapes of the LDOS and EF for the 
random model, one can see that they are clearly differ- 
ent, in contrast to the dynamical model. As mentioned 
above, in Ref.a such a difference was directly connected 
to the localization of eigenstates in the energy space. As 
a result of this localization, the spectrum statistics dif- 
fers from that predicted by the RMT. In particular for 
Wigner Band Random Matrices, the level spacing dis- 
tribution was founda to deviate from the Wigner-Dyson 
dependence. 

To check these predictions, we have calculated the 
level spacing distribution for both dynamical and ran- 
dom models for the part of the spectrum corresponding 
to chaotic eigenstates. As we expected, for the dynami- 
cal model we have observed a very good correspondence 
to the Wigner-Dyson dependence. Surprisingly, we have 
found that random model gives the same result. This 
means that the level spacing distribution is quite insen- 
sitive to the small number of localized eigenstates. On 
the other hand, this result indicates that, in the case 
of realistic matrices, the degree of level repulsion is not 
so clearly affected by the difference in the shapes of the 
LDOS and EF, as it was in the case of WBRM. 



VII. SUMMARY 

In this paper we have studied a dynamical model with 
two interacting particles (rotators). The classical version 
of this model manifests both regular and chaotic motion, 
depending on the total energy: both at small and at large 
energies the motion is regular, while, at intermediate en- 
ergies, chaotic properties are very strong. The quantum 
analog of this model can be assumed to describe two in- 
teracting spins. Our choice was restricted to the subset 
of symmetric states, which corresponds to particles with 
integer spins. 

This model has already been under investigation, both 
in the classical and in the quantum description (see, for 
example Rcf.o). However, here we have used an ap- 
proach which seems much more instructive: in the quan- 
tum case, we have represented the Hamiltonian matrix 
in the basis defined by the two-body eigenstates of the 
non-interacting system, reordered according to increas- 
ing total energy. Such a representation corresponds to 
a well known procedure in atomic and nuclear physics 
( "shell-basis representation" ) , ancLaeems to be very use- 
ful in view of recent developmentsffl. 

In this representation the Hamiltonian matrix turns 
out to be banded, with many zero elements inside the 
band. If pseudo-randomness of non-zero off-diagonal ele- 
ments is assumed (in the region of classical chaos), then 
one can refer to some modern developments of Random 
Matrix Theory: in particular, to the so-called Wigner 



Band Random Matrix ensemble, which is conjectured to 
be well suited to the description of conservative systems 
with complex behavior (seen and references therein). 
However, the assumption of pseudo-randomness of ma- 
trix elements is far from obvious: checking it was in fact 
one of the major motivations of our work. 

Random matrices in the WBRM ensemble are charac- 
terized by a sharp band inside which matrix elements are 
random, independent, and identically distributed, plus 
an additional principal diagonal with increasing entries, 
corresponding to the unperturbed spectrum of the two- 
body Hamiltonian. 

Compared to WBRM, our model has two pecuharities. 
In the first place there is no free parameter of interaction 
between the particles, the only parameter which deter- 
mines the relative strength of the interaction being the 
total energy of the system. Second, our model has an 
highly degenerate unperturbed spectrum. Still, the main 
features of the model are expected to be quite generic, 
because these peculiarities are quite typical in such phys- 
ical applications as complex atoms and nuclei. 

In this paper we have|-jiH|alyzed two main issues , mo- 
tivated by recent resultsnu'El. First, we have studied the 
structure of the eigenfunctions and of the local density 
of states (LDOS) and have compared them to what is 
known for completely random models, and for WBRM 
in particular. Second, we have looked for effects of dy- 
namical localization; though we have found no significant 
evidence for such effects, our analysi|a has brought into 
light a close connection (surmised ino) between LDOS, 
eigenfunctions, and certain classical distributions, which 
can be easily found by solving the classical equations of 
motion. 

As expected, in the region of classical regular mo- 
tion, eigenstates have a regular structure themselves; 
still, classical integrability does not result in strong local- 
ization, because these eigenstates are typically quite ex- 
tended over the basis of two-particle unperturbed states. 
On the contrary, in the region of classical chaos, the 
structure of eigenstates looks very chaotic itself. Never- 
theless, the size of such chaotic eigenstates can be smaller 
than the size of the basis, though the eigenstates may be 
treated as random (ergodic) ones on the scale of their 
localization. 

The dependence of the structure of eigenstates on en- 
ergy reflects their regular or chaotic nature, as it is 
chaotic itself in the latter case; e.g., fluctuations of the 
number of principal components are stronger where clas- 
sical chaos is stronger, which is the reason why a nonneg- 
ligeable fraction of eigenstates have a size significantly 
smaller than the basis size. 

Generally speaking, the global properties of chaotic 
eigenstates are quite similar to those found for WBRM, 
with one remarkable exception. In fact, for the WBRM 
ensemble, the expansion of exact eigenstates over the un- 
pertubed ones has a structure quite similar to the one 
observed on expanding unperturbed eigenstates on the 
basis of exact eigenstates. Instead, this symmetry is bro- 
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ken in our model, apparently due to the degeneracy of 
the unpertubed spectrum: a feature which is missing in 
WBRM. 

Expansion of unpertubed eigenstates on exact ones di- 
rectly leads to LDOS. In standard random matrix models 
the latter is known to be of the Breit-Wigner (BW) type, 
with the half-width given by the Fermi Golden Rule. In- 
stead, in our dynamical model the LDOS is BW-like only 
around the central peak; its tails have a much slower de- 
cay than predicted by the BW law. 

Dynamical localization effects are an extremely impor- 
tant issue when investigating the quantum mechanics of 
chaotiCrSystems. Our approach to this problem was based 
on Ref.B, where it was argued that, for the case of conser- 
vative systems, such effects are manifested by localization 
of eigenstates within the so-called energy shell, which is 
the range of energies ergodically explored by classical mo- 
tion. From this viewpoint, in order to detect localization 
(if any), one has to find the form of the classical energy 
shell, and then to compare it with the form of chaotic 
quantum eigenstates. 

Following this approach, we have defined and numer- 
ically computed classical distributions which strikingly 
correspond to the LDOS and to the average shape of 
eigenf unctions . Beside opening a new direction in the 
study of the global properties of quantum LDOS and 
EFs, this fact may be important for the quantum sta- 
tistical mechanics of isolated, chaotic systems of inter- 
acting particles, because the knowledge of the average 
shape of eigenstates gives analytical access to the distri- 
bution of occupation numbers of single particle statesQc. 
However, insofar as localization effects are concerned, 
the close agreement we have observed between quantum 
and classical distributions indicates that no such effect is 
present in our model, which is, in fact, in a deep quasi- 
classical region. Additional indications of absence of sig- 
nificant localization effects is provided by the analysis 
of the level-spacing distribution, which closely follows, 
in the strongly chaotic case, the predictions of Random 
Matrix theory. 

Finally, we have studied a Random Matrix analog of 
our dynamical model, which was constructed by leaving 
the unperturbed part of the Hamiltonian matrix unal- 
tered and by replacing all non zero off-diagonal elements 
by Gaussian random variables with the same mean and 
variance as in the dynamical model. In this way we were 
able to check to what extent quantum chaotic dynam- 
ics can be simulated by random interactions; in other 
words, we have checked the pseudo-randomness assump- 
tion. We have found that the random matrix model and 
the dynamical one are very similar in what concerns the 
global average properties of eigenstates. Nevertheless we 
have found that fluctuations of individual eigenstates are 
significantly stronger in the random model: in partic- 
ular there are many more eigenstates which are signifi- 
cantly more localized in comparison to the average size of 
chaotic eigenstates. In spite of this enhancement of the 
number of localized states, the level spacing distribution 



of the random model is still short of showing significant 
deviations from Random Matrix Theory. 

The most striking difference between the dynamical 
and the random model has been detected in the form of 
the LDOS. The LDOS of the random model drastically 
differs from that of the dynamical model, as it is quite 
close to the semicircle law, with an additional peak at 
the center. While the origin of this peak is related to a 
specific feature of our model, the occurrence of the semi- 
circle is somewhat surprising, because the general statis- 
tical properties of the random model are similar to those 
of the dynamical one. For WBRM, the semicircle law 
appears when the perturbation (that is, the variance of 
the off-diagonal elements) is strong; moreover, localiza- 
tion in the energy shell was found to appear only in the 
presence of the semicircle law. In contrast to the dynam- 
ical model, neither WBRM nor the random model have 
a classical analog (although the latter is much closer to a 
realistic systems than WBRM); therefore one can ask the 
question, whether the semicircle law for the LDOS can 
appear at all in quantum systems with a chaotic classi- 
cal limit. Our analysis shows that great care has to be 
taken in extending predictions of Random Matrix The- 
ory to systems of the latter class, at least if the systems 
themselves are in a quasiclassical regime. In that case 
the pseudo-randomness assumption obliterates dynami- 
cal correlations to which LDOS, and similar quantities, 
are quite sensitive. 

VIII. ACKNOWLEDGMENTS 

(FMI) thanks with pleasure the colleagues of the Uni- 
versity of Milan at Como for the hospitality during his 
visit when this work was done; he acknowledges the sup- 
port from the Grant of the Cariplo Foundation for Re- 
search and partial support from the INTAS Grant No. 
94-2058. 



IX. APPENDIX 1 

In this appendix we show that assuming a continuous 
distribution of off-diagonal non zero elements, average 
and variance can be estimated semiclassically and good 
agreement with numerical data is found. Let us start 
with Eq. (^ from which we have : 

*2 M M 

i=i j=i 

where, as usual = h^l{l + 1). Integrals can be easily 
evaluated and one has : 
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(16) 



In the same way 



/ fi2 \ 2 A/ M 

^ ' 1=1 j=\ 



h 

2 
6 



-I 



in such a way that 



(17) 



(18) 



The agreement between Eqs.(|l^),(|l8|) and numerical 
data is shown in Fig.nsl 
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FIG. 15. a) Off diagonal non-zero matrix elements aver- 
age as a function of classical L and I — 9 (circles), I — 19 
(squares), I = 39 (crosses); line is the semiclassical expres- 
sion {nL/8f. b) Off diagonal non-zero matrix elements vari- 
ance as a function of classical L and / = 9 (circles), I — 19 
(squares), / = 39 (crosses); line is the semiclassical expression 



X. APPENDIX 2 



ts,s + l 



V[(^-S+l)2-j2][(? + ,)2-j-2 



(20) 



The distance between two neighboring perturbed levels 
can be estimated as the difference between two neighbor- 
ing matrix elements (due to the symmetry of the matrix), 



'■s+l,s+2 



51, 

2 



(21) 



where the approximation is taken for j — s ^ 0. This 
means that the total splitting is of the order 



2,2 



AE - 2pihy2 - nn 



(22) 



Though this approximation is obtained for j — 0, 
corresponding to the biggest iJo-shell, similar behav- 
ior is expected for other shells. The expression ( [2^ ) 
has been checked numerically, see Fig.^^ where we plot 
AE = Eu — El as a function of h^P for different I and h 
(here E.^ and Ei are the energy of the upper and lower 
split level within one i?o-shell) For the comparison, in 
Fig.|l^, the relation AE = 2(1 + h^P) is also shown, 
which, in the classical limit ?i — > 0, I —^ oo, and for 
L > 1 yields the expression for the classical energy shell, 



AE = 2 E„ 



2(1 + L2). 




It is instructive to estimate the splitting of the energy 
levels within one iJo-shell due to the perturbation, using 
degenerate perturbation theory. Let us consider, for in- 
stance, a shell with Hq — 2fij > which has degeneracy 
p = I + 1 — j . Perturbed energy levels can be calculated 
by diagonalizing the matrix: 

hs,s' = (s,2j-s|y|s',2j-s') (19) 

This is a symmetric tridiagonal matrix with zero elements 
along the principal diagonal, whose elements for any s > 
2 are given by 



FIG. 16. Energy spectrum splitting AE due to the pertur- 
bation V, as a function of h^l'^ for 0.01 < h < 2, I < 40 
and L > 1 (points). Full line is the classical energy shell 
AE = 2(L^ + 1) . 
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